Covariant Vortex In Superconducting— Superfluid— Normal Fluid Mixtures 

with Stiff Equation of State 
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Abstract. The integrals of motion for a cylindrically symmetric stationary vortex are obtained in a 
covariant description of a mixture of interacting superconductors, superfluids and normal fluids. The 
relevant integrated stress-energy coefficients for the vortex with respect to a vortex-free reference 
state are calculated in the approximation of a "stiff" , i.e. least compressible, relativistic equation of 
state for the fluid mixture. As an illustration of the foregoing general results, we discuss their appli- 
cation to some of the well known examples of "real" superfluid and superconducting systems that 
are contained as special cases. These include Landau's two-fluid model, uncharged binary super- 
fluid mixtures, rotating conventional superconductors and the superfluid neutron-proton-electron 
plasma in the outer core of neutron stars. 
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I. INTRODUCTION 

The subject of investigation in the present work is the 
structure and energy of a stationary and cylindrically 
symmetric quantized vortex in an interacting multi-fluid 
mixture, which may consist of charged and uncharged 
superfluids and of normal fluids. This analysis has ini- 
tially been motivated by the superfluid mixture com- 
monly found in neutron star models, namely in the outer 
core region, where superfluid neutrons, superconducting 
protons and normal electrons are generally thought to 
coexist. However, due to the generality of the present 
approach, it is equally well applicable to superfluid and 
superconducting systems found in more common labora- 
tory contexts, some of which will be discussed briefly in 
the concluding section VIII . 

The study of superfluid mixtures has a long history, 
beginning with the pioneering work of Khalatnikov 
later followed by the analysis of Andreev and Bashkin 

, who incorporated allowance for a (nondissipative) in- 
teraction between the superfluids. This effect is called 
"entrainment" (sometimes also "drag" ) and plays a cen- 
tral role in the study of such fluid mixtures. The model 
has been further extended by Vardanian and Sedrakian 
to include charged fluids, and later a Hamiltonian 
formulation in the Newtonian framework has been de- 
veloped by Mendell and Lindblom The problem of 
vortices in such mixtures has been considered especially 
in the context of neutron stars, namely by Sedrakian and 
Shahabasian Alpar, Langer, Sauls Mendell and 
others. 

The covariant vortex solution in a single uncharged 
superfluid has been analyzed by Carter and Langlois , 
who have also considered the modifications due to the 
compressibility of the superfluid. The present work is on 
the one hand a generalization of this analysis to arbitrary 
fluid mixtures, including charged ones and their coupling 



to electromagnetic fields, but on the other hand is re- 
stricted (for technical reasons) to the case of a "stiff" 
equation of state. This "stiff" case is characterized by 
the speed(s) of sound being equal to the speed of fight, 
and is, within the limits of causality, the closest analogue 
to the common Newtonian incompressible models. Com- 
pressibility effects will be subject of future work. Finally, 
we mention the previously found result ^ for a Newto- 
nian vortex in a rotating superconductor, that the (hy- 
drodynamic) vortex energy is strictly independent of the 
rotating "normal fluid" of positively charged ions, a re- 
sult that will be found here to hold under much more 
general conditions. 

In the present work we will consider only stationary 
situations, which has two major advantages. First, it 
restricts the normal fluids to be in a state of rigid mo- 
tion, and moreover in the same state of rigid motion, 
because normal fiuids always possess some nonvanishing 
amount of viscosity and mutual friction. This even al- 
lows to describe a solid component in the present frame- 
work as a "normal fluid", because in the rigid state of 
motion the anisotropic effects of viscosity and elasticity 
become irrelevant. So we can for example conveniently 
describe a conventional laboratory superconductor as a 
superconducting-normal fluid mixture, consisting of su- 
perconducting electrons and a "normal" lattice of ions, 
as will briefly be discussed in the concluding section. The 
second and even more powerful consequence of station- 
arity is that we can use a conservative model based on 
a Lagrangian formalism that has been developed in re- 
cent years |l^Jll|] in a generally covariant language. The 
use of a generally covariant instead of simply Newtonian 
description has also been motivated initially by the per- 
spective of application to neutron stars, where relativistic 
effects inevitably come into play, but this approach turns 
out to be generally more flexible and convenient for the 
hydrodynamic description of such systems, even if rela- 



1 



tivistic effects are not important. 

The plan of this work is as follows. In Sec. |l| we intro- 
duce the relevant notions and equations of the covariant 
multi-fluid formalism on which the present analysis is 
based. In Sec. |l| we discuss the description of super- 
fluids in this framework and the topology of the vortex- 
type configurations. Sec. IV introduces what we called 



the "mongrel" representation of superfluid-normal mix- 
tures, that consists of choosing the superfluid momenta 
and the normal currents as the basic variables of the de- 
scription, and which will be particularly convenient for 
the present problem. In Sec. ^ we specify the class of 
cylindrically symmetric and stationary vortex configura- 
tions and obtain the first integrals of motion for these 
solutions. Sec. VI is devoted to the specification and the 
properties of the reference state, needed to separate the 
quantities attributed to the vortex from the fluid back- 
ground. Finally, the relevant vortex stress-energy coeffi- 
cients are integrated in Sec. VII, using the most general 
hydrodynamic modelization for the vortex core, and we 
find that the "rotation energy cancellation lemma" of 
still holds under the more general c ondit ions of the 
present work. In the concluding section VIII, we briefly 



illustrate the application of the foregoing results to some 
of the well known examples of superfluid and supercon- 
ducting systems. 



II. COVARIANT DESCRIPTION OF PERFECT 
FLUID MIXTURES 

The general class of (non-dissipative) mixtures of 
charged or neutral perfect fluids has been shown by 
Carter to be describable by an elegant covariant ac- 
tion principle. In this section we will briefly introduce the 
part of the formalism and notations that will be relevant 
to the present work. 

In the absence of electromagnetic effects, a mixture of 
perfect fluids can be described by a Lagrangian density 
A that depends only on the particle number currents , 
where late Latin indices, X, Y etc., enumerate the differ- 
ent fluid constituents. Variation of A with respect to the 
currents. 



6k 



/4 



(1) 



defines the dynamical momenta per particle ^ as the 
conjugate variables of the currents with respect to 
A. Here and in the following we use implicit summa- 
tion (except otherwise stated) over identical spacetime 
as well as constituent indices. Legendre transformation 
with respect to the currents, i.e. 



V = k 



/4 



(2) 



defines the "Hamiltonian density" P as a function of the 
dynamic momenta /i^ . This function only exists for non- 



degenerate systems, that is, if the functions /^(ny) de- 
fined in (|^) are invertible. The conjugate relations can 
then be written as 



dV 



(3) 



Furthermore, the form of these relations is constrained by 
the requirement of covariance, namely V (as well as A) 
has to be a scalar density, and can therefore only depend 
on scalars, i.e. on /i^/i^". This restricts relation (|3|) to 
be of the form 



TT-x = KxY Ai 



Ya 



(4) 



where the (necessarily symmetric) matrix Kxv is defined 
as 

„ dV 



XY 



(5) 



The condition of a non-degenerate system is equivalent 
to det K ^ 0, and so we can write the inverse relation 



UYa 



with Kyz = 5z ■ (6) 



In the case of noninteracting fluids, the Hamiltonian V 
would not depend on crossed scalars withX / Y ^ 

but only on diagonal terms (^//^"). In this case the 
matrix Kxy would be diagonal, and each current would 
be aligned with the respective momentum, similar to the 
case of a single perfect fluid, but any interaction terms 
between different fluid constituents in the Hamiltonian 
will lead to nondiagonal components of Kxy , and there- 
fore the currents will become linear combinations (in each 
point) of the respective momenta. This (nondissipative) 
effect is called "entrainment" and has first been consid- 
ered for superfluid mixtures of '^He and ^He by Andreev 
and Bashkin 

Before we come to the equations of motion, we need 
to extend our description to include the electromagnetic 
field and its coupling to charged fluids. This is done 
via the standard "minimal coupling" prescription that 
consists of defining the total Lagrangian density C as 



£ = A + rA, 



16n 



(7) 



where we are using units with c = 1 . The electric current 
j" is defined as 



(8) 



with being the charge per particle of the constituent 
X. The electromagnetic 2-form F^p is defined as the 
exterior derivative of the gauge 1-form Aq, i.e. 



Fal3 = 2V[qA/3] 



(9) 



where square brackets indicate (averaged) index antisym- 
metrization. The symbol Vq. denotes the usual covari- 
ant derivative, but we note that because of the antisym- 
metrization, exterior derivatives are independent of the 
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afRne connection, so we could as well replace Vq by the 
partial derivative da- 

The conjugate variables of the currents with respect 
to the total Lagrangian C are the canonical momenta tt^ , 
defined as 



dC 
dn 



(10) 



■X 



which can be seen from (^ and (Q) to be directly related 
to the dynamical momenta ^ , namely 



(11) 



The equations of motion are to be derived from the to- 
tal Lagrangian C via an appropriate variational principle. 
Imposing invariance of the action under free (infinitesi- 
mal) variations of the gauge field Aa leads to the Maxwell 
source equation, 



(12) 



However, the equations of motion for the fluids can not 
be derived via free variations of the currents this 
would simply lead to the trivial equations tt^ = 0. This is 
because free variations of the currents contain too many 
degrees of freedom, which results in overdetermined equa- 
tions of motion, therefore the variations have to be con- 
strained. It has been shown in [ pO{ that variations 6n^ 
with the correct number of degrees of freedom are gen- 
erated by infinitesimal displacements of the worldlines 
of fluid particles. These worldline variations satisfy the 
physical constraint of conserving the number of particles, 
and they result in the correct equations of motion for the 
fluids. Without entering into the technical details of this 
procedure (see ]lO| , pT[ |), the resulting equation of motion 
for each fluid X is found as (no sum overX) 



(13) 



and by contracting this equation with , we see that it 
implies that the currents are conserved, i.e. VqTiJ = 0, 
so the equations of motion reduce to the simple form of 
a vorticity conserving flow, namely (no sum overX) 



(14) 



where the (canonical) vorticity 2-form ^v^|^ is defined as 
the exterior derivative of the canonical momentum tt^ , 
i.e. 



2Vr 



77^ 



(15) 



The very compact form ( |l4|) of the equation of motion 
can be seen to "reduce" in the nonrclativistic limit to the 
(much less compact) Euler equation of a charged fluid in 
electromagnetic fields, and possibly subject to further po- 
tential forces. This is an example that shows the advan- 
tage and convenience of the covariant formalism, espe- 
cially for more complex applications like interacting mix- 
tures of possibly charged fluids in electromagnetic fields, 
as considered in the present analysis. 



And finally, the stress-energy tensor T"^ is found |Tl| ] 
in the form 

(16) 

which (in the absence of external forces) satisfies the 
equation of (pseudo) conservation, \/aT"^ — 0. From 
the form of the stress-energy tensor (|l^) we see that V 
plays the role of a generalized pressure, which reduces to 
the ordinary pressure in the case of a single fluid. 



III. PROPERTIES OF SUPERFLUIDS AND 
TOPOLOGY OF VORTEX SOLUTIONS 

We want to allow for some of the fluids to be super- 
fluid or superconducting, and we will denote these con- 
stituents by capital Greek indices Y, W etc. For "nor- 
mal" fluids (i.e. not superfluid or superconducting), we 
will use early Latin capital indices ^, i? etc., so a sum 
over all fluids (indexed by A, Y etc.) can be written as 
J2x — + ■ Apart from the electric charge there 
seems to be no fundamental difference between super- 
fluids and superconductors, and therefore we will in the 
following refer to them as "uncharged" and "charged su- 
perfluids" respectively. We note that the present treat- 
ment considers superfluids as a subclass of perfect fluids, 
and will therefore represent some restrictions as to the 
application to strongly anisotropic superfluid phases like 
they are found in ■^He jlj], which is governed by addi- 
tional "internal" degrees of freedom like the spin and 
angular momentum of the Cooper pairs. But at least for 
situations where these additional degrees of freedom of 
the order parameter can be considered as "frozen" and 
the dynamics mainly governed by the superfluid "phase" 
to be discussed in the following, the present approach 
should still represent an acceptable approximation. 

We distinguish the (connected) spacetime domain P-*" 
occupied by the superfluid constituent T from the subset 
of its respective "superfluid domain" C , which 
corresponds to what is sometimes called the "bulk" . In 
the superfluid domain S''~ the canonical momentum ttJ 
always obeys the constraint 



(17) 



where the "phase" (p^ is a continuously differentiable 
scalar on S''~ , that can be multi-valued, but the differ- 
ences between values in the same point are restricted to 
be integer multiples of 27r. This reminds of an angle 
variable and reflects the role of (p^ as a quantum phase 
e"^. In addition to the property of (quantized) potential 
flow (|l^), the superfluid T in its superfluid domain 
is perfectly inviscid. In that sense a superfluid is proba- 
bly the best representation of a perfect fluid in nature. 
On the other hand, outside its superfluid domain, i.e. in 
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j-fT y ^ ^Yie superfluid is not constrained to potential 
flow (|l7| ) and can also possess some viscosity like a "nor- 
mal" fluid. The property (17) implies that the canonical 
vorticity w"^^ vanishes on the whole superfluid domain 
5^, i.e. 



= 2V 



ttI=0, 



(18) 



which states that the superfluid is irrotational, and im- 
plies that the equation of motion (^4|) is automatically 
satisfied on . 

Irrotational flow is of course not restricted to superflu- 
ids, and the vortex-type conflgurations to be discussed 
later have been known long before the discovery of super- 
fluids, familiar examples are tornados or water flowing 
out the drain of the bath tub. But the multi-valuedness 
of the "phase" of a perfect fluid in a state of potential 
flow is not subject to a "quantization" condition of in- 
teger multiples of 27r, and a perfect fluid only exists as 
an idealization of a "real" fluid with some nonvanishing 
amount of viscosity, contrary to the completely inviscid 
superfluids in the superfluid domain. Furthermore there 
is an important energy gain associated with the super- 
fluid domain S''^ , the so-called "condensation energy". 
Superfluids consequently try to maximize their superfluid 
domain (and thereby to satisfy dl^)) as far as possible 
in the limits of the fluid domain T>. 

One of the most important consequences of (|l^) is that 
it allows for the topologically stable flow conflgurations 
known as vortices, which are characterized by the prop- 
erty that different values of the (multi-valued) phase (p''^ 
in the same point can be connected by closed paths F that 
lie entirely in the superfluid domain S'^ . As stated above, 
the difference can only be of the form 2ttN^ , where the 
integer N"^ is called the "winding number" . The winding 
number N'^ of a closed path F C can be written as 



1 

2TTh 



(19) 



It is evident that N''~ does not change for continuously de- 
formed paths F — > F' C S'^ , and N''~ is therefore a topo- 
logical constant for each equivalence class of closed paths 
in . A nonvanishing N'^ implies that the path F C 
can not be continuously contracted to a point, because it 
would necessarily have to cross at least one point P ^ S"^ 
where the phase ip^ is not deflned, and therefore iS"'' is 
necessarily multiply connected if there are nonvanishing 
winding numbers N'^ . 



IV. THE "MONGREL" REPRESENTATION OF 
SUPERFLUID-NORMAL MIXTURES 

In the previous section we have seen that a superfluid 
on its superfluid domain is generally characterized by a 
constraint (^) on the (canonical) superfluid momentum, 
while "norma " fluids are generally more easily described 



in terms of their particle number currents. For this rea- 
son it will turn out to be extremely convenient to pass 
from the "pure" type of representation used in ^ , which 
expresses all the currents in terms of all the momenta (or 
vice- versa), to a "mongrel" representation where the su- 
perfluid currents and normal momenta are expressed in 
terms of the superfluid momenta and normal currents. 
This type of representation has for example been used 
tacitly as the base of Landau's two-fluid model for su- 
perfluid ^He ||l^, which was formulated in terms of a 
"superfluid velocity" , representing in fact the irrotational 
superfluid momentum of ( [1^ ) (divided by a flxed mass), 
and of a "normal fluid" velocity, which represents the 
real mean velocity of the viscous gas of excitations. This 
will be seen in some more detail in the discussion of the 



two-fluid model in the concluding section VIII. 

In order to pass to this mongrel representation, we 
decompose the entrainment matrix Kxy into a purely 
superfluid symmetric matrix Sr^, a symmetric matrix 
Vab of purely normal ("viscous") fluids and a "mixed" 
superfluid-normal matrix Mxa, so (^) can be written in 
this decomposition as 



nr ~ St^ fj, + Mtb 
riA = Ma^ fi"^ + Vab fJ'^ 



(20) 
(21) 



where Mtb — Mbt- For clarity we use in this section 
bold typeset for denoting spacetime vectors and covec- 
tors, as the spacetime indices are not important here and 
can be put in any consistent way. Applying the inverse 
matrix to (pl|), we can easily rewrite these relations 
in the "mongrel" form 



rir = Srif 



where we defined the new matrices 

^AB ^ , = V^Afe* 



Mta'^'^Mb^ 



(22) 
(23) 



(24) 



In this representation it is easy to see that terms of the 
form nx , e.g. in the stress-energy tensor (|l^), can be 
written in the "quasi separated" form 



nxfjf^ = fJ.^ Sr.^ + riA ub ■ 



(25) 



where the effect of "mixed" entrainment between super- 
fluids and normal fluids is hidden in the use of the matrix 
S. As we consistently wrote lower constituent indices for 
currents and upper constituent indices for momenta, we 
can now use this convention to introduce a very conve- 
nient and suggestive notation, namely to use Srs' to lower 
superfluid indices T, 'P etc., and V"*^ to raise normal fluid 
indices A, B etc. This can formally be understood as 
choosing § and V as the metric tensors in the respec- 
tive constituent vector spaces of the superfluids and the 
normal fluids, but can also just be seen as a shorthand 
notation for 
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/xr = Sr.^ , and = V"^ . (26) 

In this notation, stress-energy contributions nx jJ^ take 
the simple and concise form 

nxiJ^ ^'nf^nA + fj-^ fir , (27) 

where all the information about entrainment has been 
encoded in the respective metrics of the superfluid and 
normal constituent spaces. 

We note that the superfluid constraint (|l^) generally 
applies to the canonical momenta tt^ , which only in the 
case of uncharged superfluids coincide with the dynam- 
ical momenta fi^ = "k^ — A. This implies a qualita- 
tive difference between charged and uncharged superflu- 
ids, and it will be useful to separate the superfluid con- 
stituent space into the two orthogonal subspaces that are 
naturally defined by the superfluid "charge vector" with 
components e'^ . The respective subspaces are defined 
by parallel and orthogonal projection via the projection 
tensors 

r 

Vi' = j^i^ 7 7^ = (5^ - j]^ , (28) 

where again we have used the notation er = Sr^e'^ . Now 
we can decompose constituent vectors, e.g. the superfluid 
momenta as fi^ ~ fij^ + fi^ , where 

fij = , and /x^ EE . (29) 

The subtlety of this notation is that even though a "par- 
allel" constituent vector fij^ only has nonvanishing com- 
ponents for charged superfluid constituents, and respec- 
tively fj/^ only for uncharged superfluids, the values of 
the respective components may depend on all the other 
superfluids and normal fluids, as the projection tensors 
contain the entrainment matrix S. 



V. THE STATIONARY CYLINDRICAL VORTEX 
CONFIGURATION 

In this work we will consider the simplest, because 
maximally symmetric type of vortex configuration, which 
is characterized by both stationarity and cylindrical sym- 
metry. This means that there are three independent, 
commuting (in the sense of Lie brackets) symmetry gen- 
erators fc", and m", which can be taken to corre- 
spond to time translations, longitudinal space transla- 
tions (along the vortex axis) and axial rotations, respec- 
tively. The geometric picture of the symmetry surfaces 
generated by A;", P and m" are cylindrical hypersur- 
faces that build a well behaved foliation of spacetime, 
and can therefore be parametrized by a "radial" coordi- 
nate r. Let us introduce the corresponding cylindrical co- 
ordinates ^} = {t,z,(p,r}, adapted to these 
symmetries, i.e. 



fc" = {1,0,0,0}, r = {0, 1,0,0}, m" = {0,0,1,0}. 

(30) 

The symmetry requirements and the property of con- 
served currents (p^), i.e. Vo-n^ — 0, restrict the flow to be 
purely helical, i.e., to have no radial components. There- 
fore the currents are confined to timelike hypersurfaces 
generated by the symmetry vectors and can be written 
as 

n^={4(r), nl(r),n|(r),0}. (31) 

A further consequence of the symmetry is that any phys- 
ically well defined quantity Q of the flow must be invari- 
ant under symmetry translations, which means that the 
corresponding Lie derivatives must vanish, i.e. C^Q = 0, 
for ^" being any linear combination (with constant coef- 
ficients) of the symmetry vectors fc", m" and I". This 
also holds for gauge dependent quantities like the canon- 
ical momentum tt^, provided we fix the gauge in a way 
that respects the same symmetries, i.e. when {C^A)a = 0. 
Such a gauge choice is given by 

A^^{At{r), A^r), A^ir),0}. (32) 

The components At and A^ are still subject to the resid- 
ual gauge freedom of an additive constant, i.e. 

At^At + Gt, A,^A, + g,, (33) 

but because ip is an angle variable, corresponding to 
a compact dimension, the gauge of the axial compo- 
nent Aip is completely fixed by (p2|). This is most 
easily seen by applying Stoke's theorem to a {r, (p}- 
surface integral over Fa0, i.e. J di:°'f^Faf3 = § dl'^Aa, 
which in this trivial symmetric case just reduces to 
/p °° dr (dA^/dr) = Aip{r oo), and so the gauge is fixed as 

A^(0)=0. (34) 

With the gauge choice (|3^), the symmetry condition for 
TT^ reads 

(£e^)^ = 0, (35) 

where f " can be any linear combination of the three sym- 
metry generators. The well known Cartan formula for the 
Lie derivative of a p-form Waf}-y..., namely 

{C^w)al3y... = (P + 1) 6^V[AWa/3-y...] +pV[a i^^Wxfij,,,]) , 

(36) 

can be applied to the 1-form tt^ in (^), and so we obtain 
the explicit symmetry condition. 

For superfluids (in the superfluid domain), the first term 
vanishes because of the irrotationality property (^ , and 
so the second term provides us with three independent 
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integrals of motion, corresponding to the three symmetry 
generators, namely 



- E 



r _ 



L — I TT„ 



(38) 



interpretable respectively as the energy, (canonical) lon- 
gitudinal momentum, and (canonical) angular momen- 
tum per particle. While E'^ and L"^ are generally sub- 
ject to the residual gauge freedom ( p3| ) of an additive 
constant (except in the uncharged cases e"*" = 0), the ax- 
ial constant M"^ is not, because there is no gauge freedom 
for A^p . In order to calculate the winding numbers N'^ of 
the vortex by (p^, we have to choose a path F enclosing 
the vortex axis. Such a path can always be continuously 
deformed into a path generated by m" alone, and so by 
(pq) the integration simply yields 



N' = 



(39) 



Therefore the constant (canonical) angular momentum 
per particle, M"*", is an integer multiple of H, the fun- 
damental quantum of angular momentum, and the cor- 
responding angular momentum "quantum number" is 
just the winding number N''~ . The superfluid canonical 
momenta tt"*" = ijT -\- e^ A are thereby completely deter- 
mined (in the superfluid domain) by the integrals of mo- 
tion (p8) (modulo the gauge freedom (p3)), namely 



.1 = {-E\L^ 



,nA^^,0}, with 



iv^ e z. 



(40) 



where the vanishing of the radial component tt^ follows 
from the helical direction (|3l| ) of the currents n^, and 
the entrainment relation (^ together with ( pi] ) and the 
gauge choice (^. 

In a more realistic treatment, the normal fluids are ex- 
pected to have some amount of viscosity, in which case 
the condition of stationarity, which excludes all dissipa- 
tive motion, restricts all the normal currents to be co- 
moving with the same uniform rotation 17, i.e. 



with 



A" -^f^m" = {1,0,17,0}. 



(41) 



We could also have allowed for a constant longitudinal 
velocity along P, but this is trivially annihilated by a 
Lorentz boost, and so we have chosen our reference frame 
at rest with respect to the longitudinal motion of the 
normal fluids. The symmetry condition (^5|) along the 
flowlincs of the normal fluids, i.e. with ^" = z;", together 
with the equation of motion ( p^ yields one integral of 
motion for each normal fluid, namely 



E^ = . 



(42) 



With the given restrictions on the currents (|31|) and (41) 
the integrals of motion E'^ , , N"^ , E^ and $7 are suf 
ficient for the equations of motion (Oh to be satisfied 



But in order to actually integrate these differential equa- 
tions, one is still left with the generally nontrivial prob- 
lem of solving equations for the spacetime metric gap^ 
together with Maxwell's equation ( p2|) for the gauge field 
Ac,. However, for most vortex applications of practical 
interest (including those in neutron stars), the gravita- 
tional self-interaction of the vortex can be completely 
neglected, so the background metric can in any case be 
considered as given in advance. Furthermore, as the ra- 
dial dimensions of vortices are generally much smaller 
than the lengthscale of gravitational curvature, the local 
spacetime metric of the vortex can safely be considered 
as flat, and so in cylindrical coordinates we can write it 
as 



9al3 



dx" dx'^ = -dt^ + dz^ + dip^ + dr^ . (43) 



The remaining differential equation to be solved is (|lz 
for the electromagnetic gauge field Aa. The necessary 
coefficients of the metric connection can easily be calcu- 
lated for the flat metric (p3|), and we find the explicit 
Maxwell equations for the gauge field Aa in the form 



(rA't)' = Anrj' , - {rA'j' = Anrf 



A' 

-£- \ = Anrj-^ 
r 



(44) 



(45) 



where the prime denotes differentiation with respect to r. 
Equations (^J) describe a radial electric field A[ created 
by the charge distribution j*, and an axial magnetic field 
A'^ around a longitudinal current j^. These equations 
will result in exponentially "screened" solutions, typical 
of charged superfluids. As we saw in section [II, the 
vortex is characterized by nonvanishing winding numbers 
N^, which by (0) and (|) are seen to be directly related 
to the axial components and will result in a screened 
longitudinal magnetic field Bz, which is conventionally 
defined as 



4' 

r 



(46) 



VI. REFERENCE STATE AND VORTEX 
PROPERTIES 

A. The reference state 

In the previous section we have completely specified 
the fluid configuration containing a vortex, but in order 
to separate the quantities attributed to the vortex from 
the fluid "background", we first have to specify this ref- 
erence "background" state, which will be denoted by the 
subscript 0. For any quantity Q, the part 6qQ attributed 
to the vortex is defined as the difference with respect to 
the corresponding reference value Qq, i.e. 
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(47) 



The reference state should respect at least the same sym- 
metries as the vortex state, and can therefore, by the 
reasoning in Sec. M, be characterized completely by con- 
stants i?Q, Lq, N^, Eq and Qq. Furthermore, we nat- 
urally want the the reference background to be "vortex 
free", which means that the topological constants char- 
acterizing a vortex have to vanish, i.e. Nq = 0. Another 
natural prescription is that the uniform rotation of the 
normal fluids should be the same in the reference state 
as in the vortex state, i.e. = f2. However, there is no 
such "natural" choice for the remaining constants Eq, 
Lq and Eq , if one allows for compressibility of the flu- 
ids. The compressibility is described by the fact that 
the entrainment matrix (|^) is in general a function of 
the momentum scalars, i.e. Kxy — Kxy{ij^IJ^°'), and 
therefore, if (/iJ^/i^")Q ^ 11^11^°', this generally entails 
that Kxy 7^ ^xy ■ Now, if we consider for example the t 
component of the relation (^ between currents and mo- 
menta, and if for illustration we suppose for a moment 
that there are no normal fluids, then 12^ ~ Kr^fi^* , and 



're 



K^^^Q. Choosing for example the straightfor- 
ward reference constants Eq — E^ and Lq = , leads 
to changed particle densities n^-g ^ n^, and especially 
changed mean particle number densities (in the region 
of integration with the upper cutoff radius r^c), i-c. 
rij-Q ^ rij-- We see that with this choice of reference 
constants, we compare a vortex state with a reference 
state that does not have the same number of particles in 
the region of integration. Another physically interesting 
choice of reference state would therefore rather consist in 
readjusting the reference constants Eq in such a way as 
to obtain the same mean particle number densities (and 
therefore total number of particles in the region of in- 
tegration) in the reference state. These different choices 
have been analyzed and properly accounted for in for 
the case of a vortex in an uncharged superfluid, and are 
found to be inequivalent to each other, even in the limit 

Too 00. 

Due to the additional complications of multiple en- 
trainment and charged fluids in the present analysis, we 
will postpone this problem of compressibility effects to fu- 
ture work, and restrict our attention here to the simpler 
case of a "stiff" equation of state that is characterized by 
a constant entrainment matrix, i.e. 



dK- 



XY 



= 



K^Y = Kxy 



(48) 



In this "stiff" case, the most natural reference state is 
unambiguously characterized just by choosing the longi- 
tudinal superfluid momentum components -Eq , Lq to be 
the same as in the vortex state, i.e. 

^ps{-£;^,L^,0,0}, (49) 

while the constants E^ can be fixed by taking the normal 
particle densities to be unchanged with respect to the 
vortex state, i.e. 



'Me = 'M 



n>", where = {1, 0, 0} . (50) 



Due to the assumption of a stiff equation of state ([4§|), all 
longitudinal current components rt^ and remain un- 
changed in the reference state. Furthermore we will as- 
sume the electric current to vanish in the reference state. 



0, 



(51) 



which implies that the longitudinal electric current also 
vanishes in the vortex state, 

j*=j^=0, and f=fQ=Q, (52) 

and so we also have from (|4j) (in an appropriate gauge) 



At = = , and 



0. 



(53) 



The reference state is now completely fixed by the prop- 
erties (^, ( pO| ) and (|5l|). The vortex modifies only the ip 
components of currents and momenta, so it will be con- 
venient to introduce for covectors Qa the short notation 
Q = 5qQip for the part of the that is due to the vor- 
tex, and Qq = for the part that is still present in the 
reference state, e.g. 

fil = Jl^ + fil , and A^ = A + Aq. (54) 



From (^) and (49) it is easy to see that 



A , and /Iq = —e A 



(55) 



B. The London field 

Contrary to the longitudinal components Af and Af , 
the axial gauge field Aq in the reference state will not be 
trivial, due to the uniform rotation of the charged normal 
fluids. The Maxwell equation ( ^ ) for the (p component 
in the reference state, i.e. {Aq/t)' — 0, allows for a uni- 
form magnetic field Bq in z direction (defined as in (^6|)), 
namely by integration and using (jsj) one gets, 



A' 2 
Bq = — = —Aq ^ const. , 



(56) 



where Bq is in fact the well known uniform London field 
of rotating superconductors. An explicit expression for 
the London gauge field Aq can be obtained simply from 



1^ - 
xe 



the reference property Jq = 
the "mongrel" entrainment expression (^ 
which yields 



0^ together with 
), and relation 



A, 



4) 



•'■A ■> 



(57) 



and after using (pi] ) to write = , we get the Lon- 
don field Bq as 
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Be = 2f)(e'^e^)-i (e^ 



4) 



(58) 



The London field Bq is seen to be proportional to the 
uniform rotation of the normal fluids. If we now use 
the additional property of the vanishing charge density 
( pT| ) in the reference state, i.e. Jq = 0, then we can finally 
obtain the very simple expression for the London field, 



Br 



-2n 



-2n 



erE^ 



(59) 



where we have used the notation of lowering and rais- 
ing constituent indices via the matrix § introduced in 
Sec. If we consider in particular the case of a single 
charged superfluid with mass per particle m and charge 
per particle e, this expression in the Newtonian limit. 



where E 



reduces to the well known expression 



Bq — —2VLm/e. The question of whether m in this for- 
mula should represent the bare mass or some "effective" 
mass per p articl e will be discussed briefly in the conclud- 
ing section VIII , 



C. The Magnetic field of the vortex 



The reference state properties ( }49| ) and (|5C|) further 
allow us to rewrite the axial current j"^ in the form 



r = e 
pact form 



and with (123) we obtain the com- 



— e hrir /I 



(60) 



Inserting this into the corresponding Maxwell equation 
(H) gives 



5j ~<i 



An 



-B'. 



(61) 



which can be written more explicitly as a differential 
equation for A, containing the winding numbers N'^ as 
parameters, namely 



{e'^e^)A = h{e^N'^) 



in 



-B'. 



(62) 



where the longitudinal magnetic field of the vortex, 
B = SqB^, is defined following (|6|) as B{r)j^ A'{r)/r. 
This second order differential equation for A (or B) is 
of the modified Bessel type, and the asymptotic behav- 
ior of the solutions in the limit r —> oo can be derived 
directly from this equation, namely (where means 
asymptotically proportional) 



B ^ B' 



Van A 

r — *oo 



T 



(63) 



where i is the so-called London penetration depth, which 
is given by the expression 



Aire . 



(64) 



In the Newtonian limit of a single superfluid with charge 
per particle e, mass per particle m and a particle number 
density n, the matrix § reduces to n/m, and (|6j) reduces 
to the standard expression = Aire'^n/m. 

The total electromagnetic flux of the vortex, 
<i> = / Aadx", for a circuit at sufficiently large radial dis- 
tance, is easily seen from (pSf) to be given as 



$ = 2TTh 



(65) 



which again reduces to the standard expression 
$ = N{2Tih/e) in the Newtonian limit of a single charged 
superfluid with charge per particle e. The explicit solu- 
tion of equation ( p2| ) is expressible in terms of the (mod- 
ified) Bessel functions Kq and ifi, namely 



B{r) = CaKo{rll) 



(66) 



Mr) 



2tx 



CarlK^{r/l). 



This solution is only valid in the "common superfluid 
domain", i.e. in {^j.S^ , where all the constant winding 
numbers are deflned. From the divergence of B{r) on 
the axis it is evident that the common superfluid domain 
must have a finite separation, ^ say, from the axis, which 
can be used to define what is usually called the "vortex 
core" , with ^ being the "core radius" . The constant of 
integration Cq is to be determined from the matching of 
( |66| ) with the "inner" vortex solution, i.e. for r < ^. By 
integrating (|6^) for r > ^, we get the vortex flux outside 
the core, i.e. $ — $corc, and so Co can be expressed in 
terms of the quantities ^ and the core magnetic flux ^coro, 
namely 



Co 



2T:i'^xaKi{xQ) 



(67) 



where xq is the rescaled core radius, xq = ^/^, which cor- 
responds to the inverse of the Ginzburg-Landau param- 
eter K = ll^ of the Ginzburg-Landau model. The limit 
of an extreme type-II superconductor is characterized by 
K — > cx), i.e. xq 0, Xi^K\{x{)) — > 1, so the core structure 
becomes negligible, $corc ^ and we get 



Co 



2^ 



for 



(68) 



VII. THE VORTEX ENERGY 

In this section we will consider the "macroscopic" prop- 
erties of the vortex, namely its total energy per unit 
length and the tension of the vortex line. These quanti- 
ties are obtained by integrating the local stress-energy 
tensor of the vortex, SqT^^, over the spatial section 
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{r, if} orthogonal to the ( "longitudinal" ) vortex symme- 
try axes, whose coordinates are the subset {x^} = {t, z}, 
for {i} = {0, 1}. The local stress-energy coefficients of 
the vortex are seen from (nm to have the form 



+ 



<5e (F^^Ff,^) 



g%. (69) 



The "sectional" {r, (/sj-integral is only meaningful for 
quantities that are scalars with respect to the sectional 
coordinates r and ip, and so we have to consider only 
the "longitudinally" projected tensor SqT^j. Another 
"sectional" scalar of the stress-energy tensor is the trace 
of the orthogonally projected components, which defines 
the local lateral pressure 11 of the vortex. 



211: 



rj-il 



(70) 



In the case of a "stiff" equation of state (p8|), the Tay- 
lor expansion of Vip^yX"') around the reference state 
value Vq = 'P((^^^")e) has only two terms (using (||)), 
namely 



Ya 



)• 



(71) 



The mongrel representation (Sec. |^) is particularly con- 
venient to evaluate contributions of this type, because by 
the reference property ( po|) we have i5q (n^V^n^) = 0, 
and so we find, using (^) and ([25|), 

KxY5Q{^l''^l^) = 5e{nx^l'') = Sr^JelM'^M*) • (72) 



The relevant contributions (|69|) of SqT"/^ are now 
straightforward to evaluate, and are found to be given 
by 



0, 



0, 



(73) 



(74) 



(75) 



5eV = -^Sr^ {V^ + S^'' • (76) 

Putting these results into the expression for the vortex 
stress-energy tensor (|^ , we find that the longitudinally 
projected tensor 5qT^ j is proportional to the unit tensor, 
i.e 



with 



1 

l"67r' 



T = —Se {F^^F^p) - SqV , 



(77) 



(78) 



and so the vortex energy density, SqT^°, is equal to 
the (local) longitudinal tension of the vortex, —SqT^^, 
a property that is characteristic of the stiff equation of 
state (|4^). The vortex energy per unit length U is defined 
as the sectional integral 



U 



-2tt 



dr r SqT^o = 2tt drrT 



(79) 



The energy density T can be decomposed into two parts, 

T = rvort+rrot, (80) 

where Tvort is the part that is independent of the rotation 
Q. of the normal fluids. 



1 



1 



Tvort — TT^MfM + , 



(81) 



while Trot is proportional to Vl (via Bq, see equ. (p9|)), 
Trot = Bq (^B - \e^7^A , (82) 



and the lateral pressure H, defined in ([70[), is found to be 
given by 



87r 



B^ + 2BBr 



(83) 



Expression (^Tj) for Tvoit can be transformed using 
Maxwell's equation ( |6l| ) into the "nearly integrated" 
form 



2r2 



erN^ h 



■B' 



87rr V ' 



(84) 



The easiest way to see this is to first expand only one 
in ( pl| ) using ( [55| ) and apply (61), then expand the 
remaining J? and use the second form of Maxwell's equa- 
tion (p2[) for A. In order to regroup the derivatives, one 



also has to expand one B as A! /r in the last term of (|8 
In a similar way, Tj-ot can be reduced to 



Trnt — 



^(r^B)' 



(85) 



As anticipated from the divergence of the magnetic field 
( |66| ) on the vortex axis, we encounter the same problem 
in the energy density (p4|). This well known fact is due 
to the constant superfluid (canonical) angular momen- 
tum per particle, tt^ = HN'^ , in the superfluid domain 
. Therefore each superfluid with a nonvanishing wind- 
ing number N'^ ^ 0, must have some finite "core" region 
separating the respective superfluid domain from the 
vortex axis. The actual size of the respective core re- 
gion is determined by a trade-off between the loss of 
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condensation energy associated with the core region, and 
the diverging energy density in the superfluid do- 
main. The detailed description of this superfluid-normal 
transition would ask for either a microscopic theory, or 
at least some phenomenological, e.g. Ginzburg-Landau 
type description of the involved superfluids. However, 
such detailed descriptions turn out to be unnecessary for 
our present purpose, as we can proceed on the basis of 
a very general hydrodynamic description of the vortex 
core, based only on the necessary "minimal assumptions" 
needed to avoid the energy divergence. Namely, as the 
superfluid constraint ( [l7| ) does no longer apply in the re- 
spective "core" regions, the (canonical) angular momen- 
tum ttJ there is not quantized, and is allowed to depend 
on the radial variable r. The winding number N'^ is 
strictly speaking not defined in the core region, but we 
can keep the same symbol as a shorthand notation for 
TT^/ft, so we cast our general description of the core re- 
gion in the simple form 



AA^(r) 



for 
for 



r > ^ 



(86) 



where Af^ (r) is a continuous, monotonic function, which 
has to ensure the vortex energy density T to remain finite 
on the vortex axis, i.e. in the limit r — > 0. Note that the 
"core radius" ^ is defined, as in Sec. as the radial 
distance of the "common superfluid domain" p|y. S''~ for 
the vortex axis, and is therefore the maximum core radius 
of the individual superfluids. This obviously does not 
restrict the generality of the core description (^ , as the 
Af^ (r) are allowed to remain constant until some smaller 
radius < ^. In order to have a regular behavior of the 
energy density T near the axis, it is sufficient to demand 
that A^(r) and A{r) vanish on the vortex axis at least as 



and A 



for r ^ , 



(87) 



where by "~" we mean "asymptotically proportional" 
(and not necessarily equal). This phenomenological de- 
scription is based on only two parameters, the "core ra- 
dius" ^ and the core condensation energy per unit length 
C^con- These two phenomenological parameters would 
have to be determined either from experiment or from 
a microscopic theory, but the model is now sufficiently 
determined to allow the integration of the vortex energy, 
without the need of further assumptions concerning the 
underlying physical processes of superfluidity. 
The total vortex energy per unit length is 



U — Ucon + CA'ort + C^rot , 

where according to (^9|) and (^) we have defined 

t/vort = 27r / drr TVort , 
Jo 

Urot = 27r / drr f^ot ■ 
Jo 



(89) 



The energy contribution f/iot, which is proportional to 
the rotation fl of the normal fluids, is found from (|8^) to 
be 



Urnt — 



= 0, 



(90) 



where the vanishing of the integral follows from the 
asymptotic properties ( p3| ) and (87) of the magnetic field 
B. In the Newtonian description of a rotating supercon- 
ductor, the vortex energy was already found ^1 to be 
unchanged by the rotating charged background, and this 
lemma is seen here to still hold under quite general con- 
ditions: 

Rotation energy cancellation lemma: The "hydro- 
dynamic" energy per unit length (i.e. excluding the core 
condensation energy Ucon) of a cylindrically symmetric 
and stationary vortex in a "stiff" mixture of interacting 
superfluids, superconductors and normal fluids is in- 
dependent of the uniform rotation rate fl of the normal 
fluids, despite the fact that the radial distribution of the 
hydrodynamic energy density is modified by f2, as seen in 



The vortex energy contribution J/vort in ( ^9| ) is found by 
integrating (p4[), which yields 



Cavort, = Trft 



with o<c,?/<e, 



r^ 'S>B{r)) 



(91) 



where we used the asymptotic properties (£^), (pTj), and 
the (first) mean value theorem of integration with the 
intermediate values C and rj, after a partial integration in 
the core region. We recognize two qualitatively different 
energy contributions; the first one from a "global" vor- 
tex, diverging logarithmically with the upper cutoff ra- 
dius Too , which is characteristic for vortices in uncharged 
superfluids, and the second one from a "local" vortex, 
whose energy contribution has the standard "axis field" 
form $i?(0)/87r, which is typical for vortices in charged 
superfluids. 

Using the decomposition into charged and uncharged 
superfluid subspaces via the charge projection tensors de- 
fined in (|2^), we can rewrite the first term in brackets in 
the form 



(e^e^) 



(92) 



Concerning the second term in (|9l|) , if the magnetic field 
B(r) is slowly varying inside the vortex core, then we can 
approximately replace B{ri) « Bff), and use the explicit 
expression (pa) with and (Es) to write 



5(0 = Anhie^Nf) 



1 



xoKi{xo) 



(93) 
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where xq = ^/£. In the extreme type-II hmit, where 
the core structure becomes neghgible, i.e. in the hmit 
K = l/xo > 1, where 3>cor£< Kq(£JI) « ln(£/^), and 
xqKi{xq) Ki 1, equation ( |9l|) with (|93| ) finally gives the 
simple expression for the vortex energy 



Cavort = TT^l'^ Sr*' 



(iVlArf)ln^ + (Ar|fiV|f)ln| 



(94) 



This "quasi separated" form clearly shows the respective 
contributions from a global vortex and a local vortex, 
but as mentioned above, even for vortices which have 
nonvanishing winding numbers only in either charged or 
uncharged constituents, there will generally be contribu- 
tions from both terms, due to the entrainment matrix § 
involved in the projections. 



VIII. DISCUSSION OF SOME APPLICATIONS 

In order to illustrate the foregoing general results, we 
will in this section discuss some applications to well 
known standard examples of "realistic" superfluid sys- 
tems, ordered by increasing complexity. 



A. Single uncharged superfluid 

Probably the simplest case are single, uncharged 
(isotropic) superfluids like ^He. We note that vortices 
in '^He show a much richer structure than in *He (e.g. 
see |l^), due to the anisotropic type of the microscopic 
Cooper pairing responsible for the superfluidity of "^He. 
But the present approach should still be a good ap- 
proximation at least for the '^He-B superfluid |15|, be- 
cause sufficiently far from the vortex core the additional 
(anisotropic) degrees of freedom of the order parameter 
are "frozen" and the dynamics is again mainly governed 
by the phase . 



a) at T — 0: In the case of a single superfluid at zero 
temperature, the "entrainment matrix" Kxy of (^) re- 
duces to K = nP / fjP, where n" is the particle current 
and fia the momentum per particle of the superfluid. 
There are no normal fluids, so § of ( p^ ) is given triv- 
ially hy E> = K. The charge vector vanishes, = 0, and 
the charge projection tensors are trivial, so iVj = N and 
— 0. The vortex energy (^ij) then simply reduces to 



b) at T 0: In the case of a finite temperature, the sys- 
tem can be described as an effective superfluid-normal 
fluid mixture, where the normal fluid consists of the vis- 
cous gas of excitations in the superfluid. The superfluid 
and normal particle currents are rig and tIn, and their 
respective momenta per particle /x^ and fi^ , say. There 
are no charged fluids, so = and — N. The 
entrainment matrix (|^) reads 



(96) 



and is decomposed in the mongrel representation of 
Sec. [V as V = l/i^NN, and § = Kss — K^^/K^^^, so the 



vortex energy would simply be given by inserting this ex- 
pression for § into equation (91). However, in order to 
compare this result to the usual expression for the vortex 
energy in superfluids at T ^ 0, we have to link the present 
entrainment formalism to the more common language of 
Landau's two-fluid model |Q that is expressed in terms 
of a "superfluid density" ps and a "normal density" p^. 
This "translation" has been achieved in a rigorous and 
extensive manner by Carter and Khalatnikov Jl7[ , but 
for the present purpose of an illustrative example, the fol- 
lowing very simple argument should show in a sufficiently 
convincing way how to translate between the respective 
quantities. Namely, consider the total (spatial) momen- 
tum density T"' (with z = 1,2,3) of the fluid mixture, 
for which from dl9) we have p1^= = n^ii^^ + n^/i"^'. 
Using the mongrel relations ( p2| ) and (^), this can be 
rewritten as = (/x^°S) /i^' -I- (n° V) nl, . Now we in- 
troduce the normal velocity wj, = nl^/n^, which is the 
real mean velocity of the excitations, and the superfluid 
"pseudo- velocity" vl = p^'/m^^j which is not a "real" ve- 
locity in the sense of a particle transport. In the nonrel- 
ativistic limit, where tends to the constant rest mass 



of the superfluid particles, fi rris, the irrotationality 
property of superfluids ( p^ implies V^'us' ~ 0, in other 
words "rot Vg — 0" . In these variables the total momen- 
tum density now reads 



Comparing this to the orthodox expression 



(97) 



(98) 



we can identify 



C/vort =^Vn'^ln^ 



(95) 



which is the same expression as found in Q for the single 
superfluid. In the nonrelativistic limit, where p'' — > m 
and n'^ — !■ n (where m is the rest mass of the superfluid 
particles, and n their number density), we recover the 
usual expression for the (hydrodynamic) superfluid vor- 
tex energy in the zero temperature limit (e.g. see p^ ). 



Ps = (Po) and -(<)^V. 



(99) 



This is 



P 
P 



consistent with 
Pn, namely using 



the additivity postulate 
5|) we obtain the expression 



Ps 

n^pp^ + n°p"", which effectively reduces to the total 



mass density in the Newtonian limit. In the present case 
we have — > rris for the superfluid, while jjT^ 0, as 
the normal fluid is identified with the gas of excitations, 
so the total mass density reduces to p ^ nsms- 
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In the nonrelativistic limit, expression ( |99| ) yields 
S = ps/mf, and so the equation ( |9l| ) for the vortex en- 
ergy can explicitly be written as 



C/vort ^TTh N ^ In — , 



(100) 



in agreement with the well known result in Landau's two- 
fluid model (e.g. see 0]). 



B. Two uncharged superfluids 

In the next step, let us consider a vortex in a mixture of 
two uncharged superfluids, as first considered by Andreev 
and Bashkin ||^ for a mixture of ^He and ^He. Again, at 
T = there are no normal fluids, so we have 



K 



r<i' 



K33 K34 



(101) 



The charge vector vanishes, = 0, and so Nj = and 

Nl = {N^, N^}. The expression (|9|) for the vortex en- 
ergy in this case explicitly reads 

C/vort = nh^ [{N^'fKss + {N^Kii 

+2N^N*K:ii\ In ^ . (102) 

We see that there is a purely hydrodynamic interaction 
energy due to entrainment (i.e. not related to the conden- 
sation energy in the core) from the last term in brackets, 
which is either attractive or repulsive depending on the 
sign of the entrainment coefficient K^/^. 

C. Conventional Superconductors 



When we consider cases with charged superfluids, the 
simplest example is already a two constituent system, 
because a second charged component is necessary to al- 
low for global charge neutrality. This picture applies 
for example to conventional laboratory superconductors, 
where the charged superfluid (charge e~ and particle den- 
sity n_) consists of Cooper paired conduction electrons, 
while the second component is the "normal" background 
of positively charged ions (charge e"*" and particle density 
n+). In the maximally symmetric and stationary situa- 
tions considered in the present work, "normal" compo- 
nents are naturally restricted to uniform rotation (41), 



and therefore it makes no difference whether the normal 
component is actually a real "fluid" or a solid lattice like 
in the present example. 

Because of the Cooper pairing mechanism, the funda- 
mental superfluid charge carriers have to be considered 
as electron pairs, and therefore the charge per super- 
fluid particle e~ should be twice the electron charge, i.e. 
e" — — 2e, and consequently the rest mass is = 2mc, 



where rric is the electron rest mass. The entrainment 
matrix Kxy, defined in (^) can be written as 



KxY = 



(103) 



and 
tion 



the 

of Sec 



transformation into the mongrel representa- 



IV yields § = - (X_ 



and 



Nl 



l//f4._i_. The charge vector is just 
= and = N. 



and so 



The London field: In the simple case of a vortex-free 
state, i.e. with iV = 0, there is nevertheless a nonva- 
nishing uniform London field Bq if the superconduc- 
tor is rotating (rotation rate il). Equation ( ^9|) for 
the London field immediately yields for this simple case 
Bq = —2il{E/e~), where E is the energy per superfluid 
particle, i.e. E — — /Xq . If we choose a reference frame 
with L = fj,J = 0, i.e. comoving with the superconductor 
in z direction, then E can be identified with the (relativis- 
tic) chemical potential fi^ = (— /i^/i" )^/^. In the New- 
tonian limit, where /i~ « to^, the conventional Newto- 
nian chemical potential /i^jj^j^ is related to the relativistic 
chemical potential /i~ as 



H =771 



1 



/^chcm 



(104) 



where e = Mchcm/'^ ^ ^- ^^'^ London field for a rotat- 
ing superconductor can therefore be written in the form 



B. 



-2n- 



1 



/^chcm 



(105) 



It is well known that that the "entrainment" formalism 
for interacting constituents can equivalently be expressed 
in the more conventional (albeit sometimes less conve- 
nient) language of "effective masses" 0. We see that 
in the case of two-constituent superconductors, the ef- 
fect of entrainment (i.e. effective masses) cancels out in 
the expression (105) for the London field, which there- 
fore depends quite naturally on the "bare" electron rest 
mass to charge ratio m~/e~, including a "relativistic" 
correction due to the finite chemical potential M^pm of 
the electrons. We note that this cancellation only occurs 
for systems with a single superfluid constituent, where S 
is consequently a scalar and cancels out in (^). As soon 
as there is a second (interacting) superfluid constituent 
involved, as in the following example of neutron star mat- 
ter, the London fleld does depend on the effective masses 
of the constituents. We further note that the present co- 
variant treatment is intrinsically frame-independent, and 
contrary to the analysis of , we find that Bq does not 
depend on the chemical potential of the "normal" 
component of positively charged ions. 

A very crude estimate of the relativistic correction 
term M^cm/"^~ ^O'" ^ superconductor at T = 0, tak- 
ing /i^gjjj simply to be the Fermi energy of a free electron 
gas, yields a (positive) correction of the order 10^^. This 
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is in qualitative and nearly quantitative agreement with 
precision measurements performed on a rotating Nb su- 
perconductqr |l9|. But in order to effectively compare 
expression ( |l05| ) with experimental results, a more care- 
ful estimation of fJ-^^^^^ would be necessary. 

Vortices: Now let us consider a vortex configuration, 
i.e. with N 0. We see that a similar can cella tion of the 
entrainment effect as for the London field (| 105 | ) a rises for 
the total flux of the vortex, which is seen by (|65D to give 
the usual 



$ = N^o , with $0 = , (106) 



2Trh 
e~ 



while the London penetration depth ( |64[ ) is modified 
by entrainment, namely — 47r(e^)^S. To write this 
more explicitly, we note that § can be written in the 
absence of entrainment as S^"^ = and further 

§(0) = (n_/rn^)(l - (5rci), where Src\ = McThcm/™" ^^e 
same relativistic correction facto r en countered in the ex- 
pression for the London field (105). A nonvanishing 



entrainment interaction between the constituents will 
add an additional correction term ^ontr proportional to 

the matrix element , so that S can be written as 

§ = (n_/m^)(l -|- Scntr — Srci), and so the London pene- 
tration depth reads 



= 47r(e-)^ — (1 
m 



3rclJ 



(107) 



The vortex energy is given by the "magnetic" term in 
( ^ ) alone, due to A^n = N and N±_ =0, so we recover 
the usual "axis-field" expression 



$5(0) 



(108) 



which is seen in the more explicit form ( |94| ) (for the type- 
II limit, for simplicity) to depend on the effect of entrain- 
ment, namely 

C/vort = — (1 + <5e„tr - ^rcl) In ^ , (109) 

TO ^ 

but as the total vortex energy U = Uvort + C^con also de- 
pends on the largely unknown condensation energy of the 
core, the relativistic and entrainment corrections in this 
expression seem unlikely to be of observable interest. 



D. Neutron star matter (Outer core) 

In this last example we consider the case of a (cold) 
degenerate plasma consisting of neutrons, protons and 
electrons in f3 equilibrium, as relevant for the outer core 
of neutron stars (i.e. at densities > nuclear density). In 
this case one usually assumes that there is an important 
entrainment between neutrons and protons due to their 
strong interactions, while the entrainment with electrons 



is generally supposed to be negligible. We will follow this 
assumption and denote the entrainment matrix as 



XY 





Knp 







^PP 














(110) 



The calculations of the superfluid gaps for this neutron 
star matter generally suggest (see for example |^) that 
the protons will be superconducting and the neutrons 
superfluid, while the electrons remain "normal" , so this 
system would represent a superconducting-superfluid- 
normal mixture. The matrices of the mongrel represen- 
tation of Sec. rV for this system read M = 0, V = 1/Kco 
and 



K K 
K K 

^ ^ pn ^ ^ pp 



(111) 



and we define an "entrainment coefficient" 
a = Knp/Kpp. For this system the charge vectors and 
projections are nontrivial, namely 



e^r = qKpp (a, 1) 



T 









' 1 0\ 




; 




-a 



(112) 



(113) 



where q is the charge of a proton Cooper pair, i.e. 
q = 2|e|, and we further have 



NT = TV" 



1 



The London penetration depth ( |64| ) is 
and the vortex flux (p5l) is found as 



$ = (iVP -f a7V")$o , with $0 = 



2TTh 

q 



(114) 
(115) 

(116) 

(117) 



in agreement with earlier results in the literature 
The vortex energy in the type-II limit (lol) reads 



2/' Arn\2 



7rr(7V 



{Knn - a^Kpp) In + a^Kpp In 



-nh^NPfKpp In^ 



(118) 



-2Trh^N''NT?Knp In - 



Similar to the case of a mixture of two uncharged super- 
fluids, we see that the total vortex energy consists of a 
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pure n-vortex term, a pure p-vortex term (each of which 
is modified by the entrainment), while the last term rep- 
resents an attractive or repulsive (depending on the sign 
of K^p) interaction term with respect to infinite separa- 
tion. It has been suggested that the effect of entrain- 
ment between neutrons and protons could energetically 
favor a "vortex cluster" structure (i.e. a neutron vortex 
surrounded by a dense lattice of proton vortices) with 
respect to a single neutron vortex. This question can 
strictly speaking not be addressed in the present frame- 
work of perfectly axially symmetric configurations, and 
will be subject of future investigation, but the energy of 
a single n-vo rtex {N^ — 0, TV" 0) is seen from ex- 
pression ( |118| ) to be of the same order of magnitude if 
not smaller than in the absence of entrainment (a 0), 

i.e. ui°2t^'^f^'^iN''?K^nHr^/0- Any configuration 
containing more vortices is therefore rather expected to 
have a higher energy, but the possibly attractive interac- 
tion term in (118) could lead to an effective "clustering" 
of already present vortices, namely a n-vortex that "ac- 
cretes" p-vortices until saturation. 
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